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RESUMEN 

En este artfculo se presentan varfas subrutinas que han sido incorporadas 
al sistema de computacion CEINCI-LAB, para hallar la respuesta no lineal de un 
marco piano con aisladores de triple pendulo de friction; considerando para el 
efecto que la superestructura es deformable. 

La no linealidad del aislador se ha considerado mediante un modelo 
bilineal; y la solution de la ecuacion diferencial del movimiento se realiza mediante 
el Procedimiento de Espacio de Estado. 

Palabras Claves: Analisis No lineal; Historia en el tiempo; Espacio-Estado; 
Aisladores FPT. 


ABSTRACT 

In this paper it will be presented some subroutines which have been 
incorporated to CEINCI-LAB computer system, for finding the nonlinear response 
of an isolated frame with triple friction pendulum bearings. It was considered for the 
analysis a deformable superstructure. 

The nonlinear behavior of the isolators was considered by a bilinear model 
and the solution of the equation of motion was found by a state-space approach. 

Keyword: Nonlinear analysis; time history analysis; state-space FPT 
isolators. 
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1 INTRODUCCION 


Existen dos formas de realizar el analisis no lineal de estructuras con 
aisladores sfsmicos FPT. Para la primera forma de analisis, figura 1 a, se 
considera que la superestructura se comporta como un cuerpo rfgido que solo 
aporta cargas y que el comportamiento no lineal solo se produce en el sistema de 
aislacion. (Aguiar y Morales, 2016; Aguiar y Leiva, 2017). Para la segunda forma, 
figura 1 b, se considera que la superestructura si se deforma, es decir no trabaja 
como cuerpo rfgido. En este artfculo se presenta, el marco teorico resumido para 
el analisis mediante esta segunda alternativa. 


Sistema de 
aislacion 


Superestructura 



Figura 1 Modelos de analisis no lineal de estructuras aisladas: a) Superestructura 

rfgida, b) Superestructura flexible 


2 MARCO TEORICO 


Para el analisis no lineal del portico piano de “n” pisos indicado a la 
izquierda de la figura 2, se considera un modelo de masas concentradas. Como se 
puede observar a la derecha. Este modelo considera: que las vigas son 
axialmente rfgidas y el resto de elementos son flexibles, por tanto se tiene un 
grado de libertad horizontal por piso y dos grados de libertad en cada uno de los 
nudos (desplazamiento vertical y rotacion) mostrados en la parte central; y que las 
masas son puntuales y se concentran a nivel de cada piso. 
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Figura 2 Modelo de masas concentradas de un portico piano 

El sistema de ecuaciones diferenciales que gobieman el movimiento de 
sistemas de multiples grados de libertad, como el representado en la figura 2, es el 
siguiente: 


M'q + Cq + Kq = Q (1) 

Donde: M es la matriz de masas (n x n); K es la matriz de rigidez 
dependiente del tiempo (n x n); Ces la matriz de amortiguamiento (n x n) que va 
cambiando con el tiempo; Q es el vector de cargas generalizadas; q,q,q son las 
vectores de aceleracion, velocidad y desplazamiento, respectivamente 

La solucion de la ecuacion (1) puede ser obtenida usando el Metodo 
iterativo de la (3-Newmark o cualquier otro procedimiento numerico, para resolver 
las ecuaciones incrementales del movimiento. En este documento se resolvera el 
sistema de ecuaciones diferenciales del movimiento mediante el planteamiento 
espacio-estado (Aguiar, 2007). A continuacion se indica un breve resumen del 
procedimiento: 

Si se multiplica la ecuacion (1) por la inversa de la matriz de masas, M _1 , 
por la izquierda, se tiene: 

q + M'Cq + M _1 Kq = M _1 Q (2) 

Sean: 

q-q = 0 (3) 



Reemplazando las ecuaciones (3) y (4) en (2) resulta: 


Donde: 


X = FX + r 


(5) 
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0 

-M l K 


—M _1 C 


(6) 


0 

M l Q 


(7) 


La solucion al sistema viene dada por: 


Xk + 1 = AX k + Pir k+1 + P 2 (r k+1 - r k ) (8) 

En la ecuacion (8) el subindice k, corresponde al instante de tiempo k y el 
subindice k+1 al instante de tiempo k+1 . Los valores de A, P v P 2 se indican en las 
siguientes ecuaciones: 


A = e 


AtF 


P ± =F~ 1 (A-I ) 


7-1 


F~ 1 (^-P 
\ At 

-K-'C —K~ X M 


\-a) 


0 


( 9 ) 

( 10 ) 
( 11 ) 


En la ecuacion (9), At es el incremento de tiempo con el que se halla la 
respuesta en tiempo y que coindice con el incremento de tiempo con el cual se 
obtuvo el acelerograma. Esto es lo mas facil pero se puede trabajar con dos 
incrementos de tiempo diferentes, uno para el acelerograma y otro para hallar la 
respuesta en el tiempo. 

Para hallar el vector de estado X en el instante k+1 , mediante la ecuacion 
(8) se debe conocer el valor de X en el instante k, el valor de la excitacion en los 
instantes k y k+1. 

La solucion de la ecuacion (5) es la siguiente: 

r* (12) 

X(t) = exp(F t)X 0 + I exp\F(t - r)] r(r)cf(T) 

Jo 

Al discretizar la respuesta y considerando un incremento constante de 
tiempo At se puede encontrar la respuesta en un instante t = (fc + l)At en funcion 
del valor anterior para t 0 = kAt. 

r( k + l)At 

X(k At + At) = exp(F A t)X(k At) + exp{[(k + 1)A t - r ]F] [r(T)]d(r) 

^ k At 

Se considera que la variacion de la excitacion r entre el instante de tiempo 
kAt y el instante de tiempo (fc + l)At es lineal. Luego: 
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X{k At + At) = exp{F A t)X{k At) + f^ 1)M exp{[(k + 1)A t - t]F) [r(r)]d(r) 

Al sustituir esta variation de r(r) en la integral y considerando el siguiente 
cambio de variable: 


p = (k + l)At — r -> dp = —dr 


Se tiene: 


X(k At + At) = exp(F At)X(k At) 


+ 


J r At 

exp(pF) 

0 


(At - p) 


r(k At + At) - r(k At) 
At 


+ r(fc At) 


dp 


La expresion anterior de puede expresar de forma condensada como: 
X(k + 1) = AX(k) + P x r(fc + 1) + P 2 [r(k + 1) - r(fc)] 


Donde como: A,P 1 ,P 2 fueron definidas anteriormente. 

El algoritmo de Espacio de Estado puede ser considerado como la 
generalizacion de la integral de Duhamel para varios grados de libertad. 

El programa pse_nl2 encuentra la respuesta No lineal en el tiempo de un 
sistema de n grados de libertad ante una accion sismica definida por el 
acelerograma. 

A continuacion se indica el proceso a seguir para el analisis no lineal de un 
portico piano mediante las distintas subrutinas creadas en CEINCI-LAB : 

Paso 1 : 


Se define la geometria de la estructura y los grados de libertad. 

Subrutina: Structure_Geom_2D_v2%Define structure geometry 


Paso 2: 

Se definen las cargas gravitatorias sobre los elementos. Subrutina: 
cg_data4 


Paso 3: 

Se genera la matriz masas de acuerdo al modelo de masas concentradas 
indicado en el figura 2. La matriz de masas, M, tendra una dimension n x 
n, donde “n” es el numero de pisos o numero de grados de libertad 
horizontales. Subrutina: masas 


Paso 4: 
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Se genera la matriz de rigidez global de la estructura teniendo en cuenta 
todos los grados de libertad (horizontales, verticales y giros) y la matriz de 
del aislador FPT mediante las subrutinas: krigidez_ASCE41_21 y 
krigidez_isolator2, respectivamente. Se ensamblan las matrices 
mencionadas y se obtiene la matriz de rigidez total. 


Paso 5: 

Se condensa la matriz de rigidez total a coordenadas laterales y se calcula 
los valores propios y los modos de vibracion. Subrutina: K_Modos. 


Paso 6: 


Se genera la matriz de amortiguamiento, C, empleando el algoritmo de 
Wilson y Penzien (1972) de la siguiente forma: 


C = 



(13) 


C,= 


2 SjWrn 

M\ 


(M 0i)(0?M) 


(14) 


M\ = 0?M0j 


(15) 


Siendo 0j el modo de vibracion i; W ni es la frecuencia natural del modo i; ^ 
es el factor de amortiguamiento del modo i. 

Subrutinas: amortiguamiento2 y Calc_zeta. 


Paso 7: 

Para cada instante de tiempo, se encuentra el vector de cargas Q a partir 
de los datos del acelerograma para el cual se desea realizar el analisis 
sfsmico. Subrutina: Mcargas 


Paso 8: 

Se halla el vector de estado X k+1 que contiene los desplazamientos 
laterales y las velocidades, q a , q a respectivamente. Se halla el vector de 
desplazamientos y giros q b mediante q b = -K bb K ba q a y por ultimo el 

vector q = [q“], Subrutina: psn!2. 
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Paso 9: 


Se encuentran las fuerzas y momentos en cada uno de los elementos. 
Para ello, una vez se tienen las deformaciones en coordenadas globales p 
se obtienen las acciones en los elementos en coordenadas globales PG 
multiplicando la matriz de rigidez de miembro k (matriz en coordenadas 
globales) por las deformaciones p. 

PG^ = IcWpW 

Se obtiene los las acciones de los elementos en coordenadas locales, por 
medio de la matriz de paso,r 2 _ 3 , mediante la siguiente expresion: 


PL = T 2 _ 3 PG 


Donde: TV, = [® »] 
indicado en la figura 3. 


con R = 


cos oc 
-sin oc 


L 0 


sin oc 
cos oc 

0 


0 

0 

1. 


, siendo oc angulo 



Figura 3 Sistema de coordenadas locales y globales 


Paso 10: 

Se comparan las acciones de los elementos obtenidas, PL, con las curvas 
de histeresis de los elementos. En la figura 4 se indica el modelo 
constitutive usado para la superestructura (elementos de acero), el cual 
contempla el efecto Bauschinger y el efecto de endurecimiento por 
deformacion (Aguiar R., Mora D., Rodriguez M., (2015)). 
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Figura 4 Modelo constitutive para los elementos de acero 

En la figura 5 en linea entrecortada se indica el modelo constitutive usado 
para los aisladores triple pendulo de friccion, se trata de una modificacion 
al modelo de las 3 fases presentado por Mcvitty W.J., and Constantinou 
M.C., (2015). Mas adelante se realiza un resumen de la obtencion del 
Modelo de las 3 Fases_Modificado empleado para el analisis no lineal. 



Figura 5 Modelo constitutive para los aisladores FTP 


Paso 1 1 : 

Como se puede observar el analisis no lineal de una estructura de multiple 
grados de libertad, es mucho mas complejo y emplea mas tiempo que un 
analisis lineal pues para cada intervalo de tiempo es necesario resolver un 
problema estatico equivalente: tf(t) Ax(t) = A PL(t), lo que lo vuelve un 
proceso iterative, para el cual se repetiran los pasos del 4 al 1 0 
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3 SISTEMA DE COMPUTACION CEINCI-LAB 


A continuation se indica la forma de trabajar con los programas de 
CEINCI-LAB, para encontrar la respuesta en el tiempo no lineal de un portico 
piano con aisladores FTP. Las unidades de trabajo, tienen que ser consistentes, 
se recomienda trabajar en toneladas y metros. 


3.1. Definicion de las secciones de los elementos: vigas, columnas 
y aisladores FPT 


Como primer paso se deben definir las secciones de los elementos 
estructurales: vigas, columnas y aisladores sismicos en el caso de existir. Para la 
definicion de las secciones se debera indicar un codigo numerico para el tipo de 
material y forma de la seccion transversal, dichos codigos se indican en la tabla 1 . 


Tabla 1 Codigos de material y perfiles tipo. 



TIPO DE 
MATERIAL 

SECCION TRANSVERSAL 

Perfil 
tipo 1 

Perfil 

tipo 

H 

Perfil 

Tubular 

Rectangular 

Maciza 

rectangular 

Maciza 
tipo L 

Maciza 
tipo T 

HORMIGON 

2 

- 

- 

- 

1 

1 

1 

ACERO 

1 

1 

2 

3 

- 

- 

- 

AISLADORES 

SISMICOS 

FTP 

1.3 

1 


3.2. Vigas y columnas de acero 


Para definir una seccion de acero, sea una viga o columna, se utilizaran 
los campos indicados en la tabla 2. En la figura 6 se indica las secciones 
transversales consideradas, con su respectiva nomenclatura. 


Tabla 2 Datos para describir una seccion de una viga o columna de acero 


Nombre 
de la 
seccion 

= 

Tipo de 
material 

fy 

fu 

Ey 

Tipo de 
seccion 

d 

tw 

bs 

bi 

tfs 

tn 

A definir 

por el 
usuario 


Segun 

Tabla 1 

Li'mite 

elastico 

del 

acero 

Li'mite 

de 

fluencia 

del 

acero 

Modulo 

de 

Young 

Segun 

tabla 1 

Alto de 

la 

seccion 

Espesor 

del 

alma 

Ancho 

superior 
de las 
alas 

Ancho 

inferior 
de las 
alas 

Espesor 

superior 
de las 
alas 

Espesor 

inferior 
de las 
alas 
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bi 




Figura 6 Secciones transversales consideradas, para elementos de acero. 


3.3. Vigas y columnas de hormigon 


La informacion para definir las secciones de las vigas y columnas de los 
elementos de hormigon armado, se indican en la tabla 3 y en la figura 7. De igual 
manera se indican tres vectores de datos por cada elemento, para el nudo inicial, 
centra de luz y nudo final. 


Tabla 3 Datos para describir una seccion de una viga o columna de hormigon 


Nombre 
de la 
seccion 

= 

Tipo de 
material 

fc 

Ey 

fy 

Tipo de 
seccion 

h 

b 

bw 

tf 

A definir 

por el 
usuario 


Segun 

Tabla 1 

Resistencia 

del hormigon 

Modulo de 

Young del 
acero de 
refuerzo 

Resistencia 

del acero de 
refuerzo 

Segun 

tabla 1 

Alto de 

la 

seccion 

Ancho 

de la 
seccion 

Espesor 

del alma 

Espesor 

de la 
alas 


d’ 

Assup 

Asinf 

Stirrup 

separation 

Av 

Asl 

Recubrimiento o alto 
efectivo a compresion 

Area de acero 
superior 

Area de 
acero inferior 

Separation de 
estribos 

Area de las 
ramas de 
estribos 

Area de acero en el medio de la 
seccion en el caso de columnas 




As_sup 

Asl 

As2 

Asjnf 


Figura 7 Nomenclatura de vigas (izquierda) y columnas (derecha) de hormigon 

armado. 
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3.4. Aisladores de triple pendulo 


Los programas de CEINCI-LAB fueron desarrollados en esta fase de 
estudio unicamente para aisladores sfsmicos tipo Triple Pendulo de Friccion FTP. 
El aislador sfsmico vendra definido por la informacion indicada en la tabla 4, donde 
se ha indicado el significado de las variables a exception de L^Lj que son las 
distancias desde el centra de gravedad del aislador hasta el punto medio de la 
viga inferior o superior, como se ilustra en la figura 8. 


Tabla 4 Informacion para los aisladores de triple pendulo de friccion. 


Nombre 
de la 
seccion 

= 

Tipo de 
material 

Tipo de 
seccion 

fy 

qy 

fu 

qu 

U 

Lj 

kv 

A definir 
por el 
usuario 


Segun 
Tabla 1 

Segun 
tabla 1 

Li'mite 

elastico 

Desplazamiento 

elastico 

Li'mite 

plastico 

Desplazamiento 
asociado al li'mite 
plastico 



Rigidez 
vertical del 
aislador 



Figura 8 Modelo de un aislador como un elemento corto, de longitud L t + Lj 


Sea k 0 la matriz de rigidez del elemento aislador y sea k 3 la matriz de 
rigidez en coordenadas globales del elemento aislador, las mismas que se hallan 
de la siguiente manera: 


k 0 


k h 0 
. 0 k v 


-1 0 l t 1 0 lj 

. 0 -10 0 10 


k 3 =rk 0 T 
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De tal manera que las distancias li,lj, sirven para hallar la matriz 
cinematica T y con ella encontrar la matriz de rigidez en coordenadas globales. El 
significado de las variables de la matriz k 0 , es el siguiente: k h es la rigidez 
horizontal; si se trabaja con las rigideces tangentes, valdra k e para el rango 
elastico y k p para el rango plastico; en cambio si se trabaja con la rigidez secante 
k h = k ef . La rigidez vertical k v se halla con la siguiente ecuacion. 


A = 


nbl 


K 


E A 

IT 


Donde: b 4 es el diametro del aislador interior del FPT; E es el modulo de 
elasticidad del acero; H es la altura total del aislador FPT. 


3.5. Asignacion de las secciones a los elementos de la estructura 


Una vez asignada la numeracion a vigas y columnas de la estructura se 
procede a asignar las secciones antes definidas a cada uno de los elementos. 
Esta asignacion se realiza en “Tablel” de la siguiente manera: 


Tablel = 


Seccion en el nudo 

inicial de elemento 1 

Seccion en el centre luz de 

elemento 1 

Seccion en el nudo final de 

elemento 1 




Seccion en el nudo 

inicial de elemento n 

Seccion en el centra luz de 

elemento n 

Seccion en el nudo final de 

elemento n 


Elemento 

1 

Elemento 

n 


3.6. Generacion de los elementos de la estructura 


Para la generacion de la matriz de rigidez y del vector de cargas 
generalizadas los datos que se deben ingresar, se resumen a continuation. 
Primero se debe indicar el numero de pisos np y el numero de nudos restringidos 
nr. Luego para general los elementos se debe crear una matriz con la information 
indicada en la tabla 5. 


Tabla 5 Datos para generar los elementos de un portico piano. 
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GEN 


N c del 

Nudo 

Nudo 

N D de 

Incre-mento en la 

Incremenla de la 

Increment de la 


elemento 

inicial 

final 

elementae a 

generar 

numeration 

numeracidn en el nudo 
inicial 

numeracidn en el nudo 
final 


Para la creacion de las coordenadas de los nudos, la informacion 
solicitada se indica en la tabla 6. 


Tabla 6 Datos para generar las coordenadas de los nudos. 


NCDE 

= 

N B del 
nudo 

Coord. X 
del nudo 
inicial 

Coord X 
del nuda 
final 

N c de nudas 
a generar 

In-crementa en la 
nu meraoion 

Incremental de la 
tongitud en X 

Incremento de la 
longitud En Y 


3.7. Asignacion de carga muerta 


Se identifica al elemento cargado y se coloca la carga uniforme distribuida, 
en un archivo, como se indica en la tabla 7. 


Tabla 7 Asignacion de carga vertical. 


Fm 


N Q del elemento cargado 


Carga distribuida existente actuando sobre el elemento 


3.8. Asignacion de carga sismica 

Se cargara el vector que contiene la aceleracion del registro para el cual 
se desea realizar el analisis sfsmico y el intervalo de tiempo para el cual se obtuvo 
el registro. 

3.9. Modelo Constitutive para los aisladores FTP 

En la figura 9, se presenta la nomenclatura de un aislador de triple pendulo 
de friccion, el cual queda definido por doce parametros geometricos y por cuatro 
parametros friccionales. En los casos convencionales: Ri=R 4 » R 2 =R 3 , 
|i2=|J3<|ii=|i4, di=d4 y d2=d3, hi=h4 y h2=h3; de esta forma las propiedades 
geometricas se reducen a seis y las friccionales a dos. Para esta condition, Me 
Vitty Y Constantinou (2015) han definido un modelo constitutive de 3 fases 
indicado en la Figura 5. 
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Figura 9 Nomenclatura del aislador de triple pendulo. 

Fuente: Fenz y Constantinou (2007) 

En el programa de analisis no lineal desarrollado, se considera el 
diagrama de histeresis 3 Fases_Modificado indicado figura 5 (en linea 
entrecortada), como que fuera bilineal. Donde el punto de fluencia F y esta definido 
por la fuerza maxima y desplazamiento maximo de la fase 1 y el punto de maxima 
fuerza F u por la maxima fuerza y desplazamiento de la fase 2. De tal manera que 
las ecuaciones que definen el modelo bilineal, son: 


q y = 2 (jUi - Hi ) Ri.eff 

W 

q + F i w 


Fy 2 R 


2 ,eff 


q u — q y + 2 
W , 

Fu = Tr ( q ~ u *) + Fi W 

z K l,eff 


Ri.eff = /?; - /i; Para I de 1 a 4. 


dl = 


Ri 


Donde q y ,F y ,q u ,F u , son los desplazamientos y fuerzas del punto de 
fluencia y del punto que se considera de maxima capacidad del aislador. Es 
importante destacar que la fase 3 no se considera en el estudio y que corresponde 
al caso de que el aislador interior llega a los topes exteriores, con un mayor 
desplazamiento lateral; de tal manera que se esta trabajando con un modelo 
conservador debido a que el desplazamiento ultimo se considera el 
desplazamiento maximo de la fase 2 y no el de la fase 3. Falta indicar que W es el 
peso total que llega al aislador; q es el desplazamiento. 

Con los parametros anteriores se definen las rigideces elastica y post 
fluencia , K e y K p , respectivamente. 
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K e 


5l 

(Jy 


K p = 


F U Fy 
Ru ~ Ry 


4 ANALISIS NO LINEAL DE ESTRUCTURAS CON 
AISLADORES FTP 


4.1. Descripcion del portico 


Se analiza el portico B y el portico E, del Bloque Estructural 1, del nuevo 
Centro de Investigaciones Cientificas de la Universidad de Fuerzas Armadas 
ESPE. En dicho portico se van a colocar sobre los pedestales de los plintos el 
aislador de triple pendulo de friccion indicado en la figura 10 y sobre los aisladores 
una losa de aislacion de las mismas caracteristicas de la losa del primer piso. Los 
coeficientes de friccion para el aislador son: j Ui = £4 = 0.07; n 2 = |U 3 = 0.01. 



Figura 10 Geometria de aislador de triple pendulo de friccion. FTP8833/1 2-1 2/8-6 

La carga muerta considerada para el analisis es D = 0.8 T/m 2 y la carga 
viva L = 0.25 T/m 2 . Para el analisis sismico se considera una carga de D + 0.25 L. 
En la figura 11 se indica la geometria de los porticos, con la carga vertical 
actuante. En la parte inferior de la figura 1 1 se indica la seccion de la columna y 
viga tipo. 




19.76 
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PORTICO B 



PORTICO E 

5.83 T/m 



SECCION TlPICA SECCION TIPICA 

PARA COLUMNAS PARA VIGAS 


Figura 11 Cargas actuantes en los porticos con aisladores sfsmicos. 
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La rigidez elastica y de pos fluencia del aislador se calculan con las 
expresiones indicadas en el la seccion 3.9. En la tabla 8 se presenta, la carga 
vertical total que gravita en cada aislador, las fuerzas y desplazamientos que 
definen el modelo bilineal indicado en la figura 5 para el portico B. 


Tabla 8 Carga vertical y valores de fuerza y desplazamiento que definen el modelo 

bilineal del aislador. 


Aislador 

W(t) 

Fy(t) 

qy (m) 

Fu (t) 

qu (m) 

1 

97.42 

6.8194 

0.0275 

18.4444174 

0.5366 

2 

182.5 

12.775 

0.0275 

34.5525168 

0.5366 

3 

86.76 

6.0732 

0.0275 

16.4261718 

0.5366 


4.2. Analisis no lineal 


Para el analisis no lineal de historia en el tiempo se considera el 
acelerograma indicado en el figura 12, correspondiente al registro de Portoviejo 
para intervalos de tiempo de 0.01s. 



Figura 12 Registro sismo de Portoviejo del terremoto del 16 de abril de 2016 

En la figura 13 los cuadros coloreados en gris, significan que los unicos 
elementos que trabajan en el rango no lineal son los aisladores mientras la 
superestructura permanece en el rango elastico. Los desplazamientos sufridos por 
sistema de aislacion durante el intervalo de tiempo considerado son indicados en 
la figura 16. 
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PORTICO B 



10 15 20 


PORTICO E 



Figura 13 Comportamiento de los porticos pianos ante la carga sfsmica 

La respuesta en el tiempo del sistema de aislacion asf como su 
comportamiento histeretico se representan en las figuras 14 y 15 respectivamente. 
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Figura 14 Respuesta en el tiempo del sistema de aislacion para los porticos 

analizados 


PORTICO B 


PORTICO E 




Figura 15 Comportamiento histeretico del sistema de aislacion de los portico 

analizados 


La respuesta para cada instante de tiempo de los grados de libertad del 
portico B y E son los mostrados en las figuras 1 6 y 1 7 respectivamente. 
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Figura 16 Respuesta en el tiempo por piso del portico B 
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Figura 17 Respuesta en el tiempo por piso del portico E 


En la Tabla 9 se recopilan los desplazamientos maximos resultantes para 
cada grado de libertad, los cuales son representados en la figura 1 8. 


Tabla 9 Desplazamientos maximos 


PORTICO B 

Piso 

Desplazamiento 
Maximo (m) 

1 

(Sistema Aislacion) 

0.31196478 

2 

0.319854166 

3 

0.326369537 

4 

0.330477331 
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PORTICO E 

Piso 

Desplazamiento 
Maximo (m) 

1 

(Sistema Aislacion) 

0.336715131 

2 

0.348615464 

3 

0.359388748 

4 

0.366530434 

5 

0.368233967 


Portico B 


Portico E 


18 

16 

14 

12 

10 

8 

6 

4 

2 

0 


0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Desplazamiento (m) 


18 

16 



0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 


Desplazamiento (m) 


Figura 18 Desplazamientos maximos en los porticos analizados 


5 Uso de CEINCI-LAB 


A modo de ejemplo del uso de los programas CEINCI-LAB, en la figura 21 
se presenta el programa, con el cual se halla la respuesta no lineal del portico B y 
en figura 19 se indica la numeration de nudos y elementos. Para los nudos se 
debe numerar primero los restringidos y para los elementos es conveniente que el 
elemento aislador se numere al final. 

Para el analisis sismico se considera que las vigas son axialmente rigidas, 
con lo que se tiene un grado de libertad horizontal por piso, los mismos que se 
numeran primero, como se observa en la figura 20; luego se numeran los 
restantes grados de libertad. Se halla la matriz de rigidez de la estructura con 
todos los grados de libertad y se condensa a las coordenadas laterales. 



Figura 19 Numeracion de nudos y elementos de Portico B con aisladores 

sfsmicos. 



Figura 20 Grados de libertad considerados en el analisis. 
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• Listado del programs 

1 % ANALISIS NO LINEAL DE HISTORIA EN EL TIEMPO DE PORTICO CON AISLADORES 

2 %TRIPLE PENDULO DE FRICCION 

3 

4 % BLOOUE 1 PORTICO B (B1PB) 

5 

6 %NOTA : TODOS LOS VECTORES DE SECCIONES DEBEN SER DE LA MISMA DIMENSION 

7 % PARA COMPLETAR LA DIMENSIONS SE PUEDE COMPLETAR CON CEROS 

8 % Unidades: toneladas y metros 

9 

10 - clear 

11- set (gcf, ' color ' , 'w‘ ) ; 

12 - Ey=2 0430000 ; %Modulo de elasticidad del acero 

13 

14 % DEFINCION DE LAS SECCIONES 

15 

16- SecC=[l, 35150, 0,Ey, 3, .55, .02, .55, .55, .02, .02] ; %columnas rectangulares tubulares de acero 

17- SecV=[l, 35150, 0,Ey,l, .55, .0125, .3, .3, .025, .025] ; %vigas tipo I de acero 

18 %Isolator: [Material (1 . 3) ,Type (1) , Fy, qy, Fu, qu, Li,Lj , lev] 

19 - Seclsl=[l. 3, 1,6. 8159, 0.0275, 18. 4350, 0.5366, .45, .45,4750000,0,0] ;%aisladores FTP 

20 - Secls2= [1.3, 1,12. 7701, 0.0275, 34. 5393, 0.5366, .45, .45,4750000,0,0] ; 

21 - Secls3=[l. 3, 1,5. 9542, 0.0275, 16. 1043, 0.5366, .45, .45,4750000,0,0]; 

22 

23 % AS I GNAC ION DE LAS SECCIONES A LOS ELEMENTOS DE LA ESTR. 

24 

25 - Tablel= [SecC SecC SecC %EXeml 

26 SecC SecC SecC %Elem2 

27 SecC SecC SecC %Eleia3 

28 SecC SecC SecC %Elem4 

29 SecC SecC SecC %EIem5 

30 SecC SecC SecC %ELem6 

31 SecC SecC SecC %ELem7 

32 SecC SecC SecC %ElemS 

33 SecC SecC SecC %Elern9 

34 SecV SecV SecV %Eleml0 

35 SecV SecV SecV %ELernll 

36 SecV SecV SecV %Eleml2 

37 SecV SecV SecV %EIeml3 

38 SecV SecV SecV %Eleml4 

39 SecV SecV SecV %Eleml5 

40 SecV SecV SecV %EIeml6 

41 SecV SecV SecV %Eleml7 

42 Seclsl Seclsl Seclsl %Eleml3 

43 Secls2 Secl32 Secls2 %Eleml9 

44 Secl33 Secls3 Secls3] ; %Elem20 

45 

46 
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47 % GENERACION DE LA GEOMETRIA DE LA E STRUCT. j 

48 - np=4;%numero de pisos 

49 - nr= [1 111 

50 2 1 1 1 

51 311 l];%numero de nudos restringidos 

52 

53 %Generacion de los elementos de la estructura 

54 - GEN= [ 1 4 7 2 3 3 3; 2 5 8 2 3 3 3; 3 6 9 2 3 3 3; 

55 10 4 5 3 2 3 3; 11 5 6 3 2 3 3; 

56 18 1 4 2 1 1 1] ; 

57 

58 %Generacion de los nudos de la estructura 

59 

60 - NODE= [ 1 0 0 1 1 10.99 0;3 20.59 0 0 0 0 0;%[NUDO, X, Y, NUDOS A CREAR, INCREME NTO # DE NUDO,dX,dY] 


61 4 0 0.9 3 3 0 4.94; 5 10.99 0.9 3 3 0 4.94; 6 20.59 0.9 3 3 0 4.94]; 

62 

63 % AS I GNAC I ON DE CARGAS 

64 - Fm= [10 4.43 

65 11 4.43 

66 12 4.43 

67 13 4.43 

68 14 4.43 

69 15 4.43 

70 16 4.43 

71 17 4.43]; 

72 % CARGAS SISMICAS 


73 — load ( ' a_POE ' ) %Registro de acelerograma 

74 - p=a_POE*9.81;%m/s2 

75 - dt=0 . 01; %intervalo de tiempo de acelerograma 

76 - plot ( (1 : length (p) ) *dt,p ’), input (' enter '), close (gcf) 

Figura 21 Programa utilizado para analisis no lineal de historia en el tiempo. 
El mismo procedimiento se realiza para el portico E. 


6 CONCLUSIONES 


En este documento se presenta el uso de las distintas subrutinas 
desarrolladas en CEINCI-LAB para la obtencion de la respuesta no lineal de 
porticos pianos aislados mediante FPT considerando la superestructura flexible. 

Como se puede observar en los resultados de los desplazamientos 
obtenidos, la diferencia entre el desplazamiento del sistema de aislacion y el tope 
de la estructura no es mas 1 .5 cm, por tanto se puede concluir que para etapas de 
prediseno la consideracion de superestructura rfgida Neva a resultados aceptables. 
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